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ABSTRACT 

Context. Lyman alpha (Ly a) emitting galaxies (LAEs) are used to probe the distant universe and are therefore important for galaxy 
evolution studies and for providing clues to the nature of the epoch of reionization. However, the exact circumstances under which 
Ly a escapes a galaxy are still not fully understood. 

Aims. The Trident project is designed to simultaneously examine Ly a, H a, and Lyman Continuum emission from galaxies at redshift 
z~-2, thus linking together these three aspects of ionising radiation in galaxies. In this paper, we outline the strategy of this project and 
examine the properties of LAEs in the GOODS North field. 

Methods. We performed a narrowband LAE survey in GOODS North using existing and two custom made filters at the Nordic 
Optical Telescope with MOSCA. We use complementary broad band archival data in the field to make a careful candidate selection 
and perform optical to near-IR SED fitting. We also estimate far-infrared luminosities by matching our candidates to detections in 
Spitzer/MIPS 24/rm and Herschel/PACS catalogs. 

Results. We find a total of 25 LAE candidates, probing mainly the bright end of the LAE luminosity function with L^ya ~ 1-15x10'*^ 
erg s“*. They display a range of masses of ~ 0.5 - 50 x 10®Mg, and average ages from a few tens of Myr to 1 Gyr when assuming a 
constant star formation history. The majority of our candidates also show signs of recent elevated star formation. Three candidates have 
counterparts in the GOODS-Herschel far-IR catalogue, with luminosities consistent with ultra-luminous infrared galaxies (ULIRGs). 
Conclusions. The wide range of parameters derived from our SED fitting, as well as part of our sample being detected as ULIRGs, 
seems to indicate that at these Ly a luminosities, LAEs do not necessarily have to be young dwarfs, and that a lack of dust is not 
required for Ly a to escape. 

Key words. Galaxies: photometry. Galaxies: structure 


1. Introduction 


Few emission lines have been studied as extensively in 
the distant universe as 12 16 A (Ly a). First suggested by 
Partridge & Peeble'sl (Il967h as a means of finding high red- 
shift galaxies, it became a signature to look for in galaxy 
surveys for decades. It was however not until deep ded¬ 
icated surveys were performed that any significant quan- 
tities of these Ly a emitting (LAE) galaxies were found 
JCowie & Hul mH). Ever since, Ly a has been success- 
fully used to find galaxies across redshi f ts z ~ 2 - 7 (e.g. 


t lOuchi et alJl2003L 12005: iGawiser et alJ 120061: lAiiki et al, 

tiGronwall et alJl2007: 

Pirzkal et alJl2007l: Finkelstein et al. 

t iNilsson et al.l 12009 

: iHaves et alJ 2010al lYuma et al. 


Ouchi et al. 120101: G uaita et a^ 20 Ht lA dams et al.l 201 lb 

Blanc et ahT 2011 : Shibuva et al. 2012 r~ Wargas et alJ f2014h . 

even probing the end of the epoch of reionization (e. 
Malhotra & Rhoads| 2004 Diikstra et alJl2007l: [Ono et al 


.lensen et alJl2013: Matthee et al.ll2014l. 


* Based on observations made with the Nordic Optical Telescope, 
operated by the Nordic Optical Telescope Scientific Association at the 
Observatorio del Roque de los Muchachos, La Palma, Spain, of the In- 
stituto de Astrofisica de Canarias. 


It quickly became apparent in the earliest studies that Ly a 
is nowhere near as bright and ubiquitous i n high redshift g alax- 
ies as was first imagined (see for example IPritched d 19941) for a 
review of non-successful Ly a surveys). Recombination theory 
tells us that roughly one third of all of the ionizing flux from hot 
stars should be reprocessed into t he Ly a emission line as the 
atomic hydrogen gas recombines dOsterbrockl[l96:^ . but some¬ 
thing is keeping Ly a from escaping most galaxies. 

It was natural to assume that interstellar dust would play 
an important role in absorbing Ly a photons, as dust absorbs 
rest-frame ultraviolet light very effectively. Indeed, early stud¬ 
ies showed a te ntative anti-correlation be tween dust content and 
Ly a emission dMeier & Terlevichlll98ll) . However, it was soon 
shown that this could not be the only factor governing Ly a 
escape, as there are examples of strongly emitting Ly a galax¬ 
ies rich in dust, and nearly dust-free galaxies where Ly a i s 
completely a bsorbe d dKunth et al.lll994ilGiavalisco et alJll996l) . 
iKunth et*^ dl998h showed the importance of the interstellar 
medium (ISM) kinematics, as the Ly a photons resonantly scat¬ 
ter on neutral Hi gas. Ly a radiation at the rest wavelength is 
optically thick already at column densities Nhi ~ lO''^ cm“^, so 
even a tiny amount of neutral hydrogen will affect the emergent 
line significantly. With multiple scatterings, the path length is 
greatly increased, and even a small amount of dust can in prin¬ 
ciple absorb a large amount of photons. However, if the neutral 
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HI gas is in rapid motion compared to the Ly a sources, for ex¬ 
ample due to supernova feedback or strong stellar winds, Ly a 
can be shifted out of resonance and escape more easily. With a 
large range of observational studies, combined with advanced 
simulations of Ly a radiative transfer through galaxies (e.g. 
Verhamme et al.l l2006i 120121: iLaursen et al.l 12013t iDuval et^ 
Mir we are only now starting to build a coherent picture of 
how Ly a is produced and escapes from galaxies. 

There are still, however, many questions that remain unan¬ 
swered. Two galaxies with essentially the same properties in 
terms of mass, star formation rate, age, gas-mass ratio, dust con¬ 
tent etc., ca n show completely different behavior when it comes 
to Ly a (e.g. lHaves et al.ll201^ . We are still not certain of which 
parameters are dominant for governing Ly a escape. As the num¬ 
ber of LAE surveys is rapidly increasing, and the volumes and 
numbers probed can now be measured in whole degrees on the 
sky and amount to thousands of galaxies, it is more important 
than ever that we truly understand the mechanisms governing 
these galaxies. 

Another outstanding question in the study of high redshift 
galaxies is the manner of escape of photons with an energy 
above 13.6 eV (<912 A), capable of ionizing neutral hydrogen. 
We refer to this as Lyman continuum (Ly C) radiation in this 
paper. As the universe was gradually reionized somewhere at 
z > 6 dFan et al.ll2006h . the natural ingredient that is assumed 
responsible are the first stars and galaxies. However, in order 
to reionize all of space, the amount of ionizing radiation escap¬ 
ing out of these galaxies must have been very large. Yet, when 
we try to measure this quantity in strongly star-forming galax¬ 
ies ("starbursts"), we see relatively modest escape fractions (see 
e.g. lVanzella et al.ll2012i and references therein). The Ly C es¬ 
cape fraction remains difficult to measure directly, as the signal 
from nearby sources is absorbed by Galactic H i gas and requires 
space-based ultraviolet observatories with high sensitivity even 
at redshifts up to z < 3. 

With these unresolved questions in mind, we started the 
project that we call Trident. In this first paper, we present the 
general outline of the project and our first results, focusing on 
Ly a emitting galaxies in GOODS-N. The H a to Ly a proper¬ 
ties and the Ly C escape will be the subject of future papers. 
In this paper, we focus on the properties of the Ly a emission, 
estimating stellar properties from spectral energy districution 
(SED) analysis. We also identify three individual candidates in 
far-infrared survey data and discuss the implications on their dust 
content. 

In section|2]we outline the general idea of Trident and how it 
is structured. In section[3]we describe the data we have collected 
and how it was reduced. Section|4]describes how we selected our 
candidates from the data. In section |5] we measure far infrared 
properties of our candidates found in the GOODS-Herschel cat¬ 
alog. In section |6] we explain how we combined our data with 
archival survey data of GOODS-N, and how this was then used 
for the SED fitting routine. In section |7] we summarize and dis¬ 
cuss the implications of these results. 

Throughout the paper, we assume Ho = 70 km s 'Mpc *, 
Qm = 0.3, Qa = 0.7 


2. The Trident survey 

We chose the name Trident because we are trying to find three 
types of emission related to atomic hydrogen from the same 
galaxies in one combined survey. The idea for a combined Ly a 
and H a study builds largely on the Double-Blind Deep project 


(iHaves et al.ll20103i^ . in which these emission lines were suc¬ 
cessfully studied in galaxies in the GOODS-South field, us¬ 
ing the FORSl and HAWK-I instruments at ESO VLT. Being 
a relatively deep study (reaching line fluxes at 5cr of 6.8 and 
7.8 X 10“'^ergs“'cm“^ in H a and Ly a, respectively), it probed 
the faint end of the Ly a luminosity function (LE) at z ~ 2 
to an unprecedented a ccuracy. A similar approach was used by 
iNakaiima et al.l (1201 2l) who studied Ly a, [On] /13727 and Ha 
galaxies at z ~ 2.2 in the SXDS field. 

Drawing on our experience gained from Double-Blind Deep, 
we started our current project. One of the drawbacks of Double- 
Blind Deep was that it was designed to find predominantly faint 
emitters. However, constraining the bright end of the Ly a LE 
requires a larger survey volume, as bright emitters are rare, and 
this is one of the goals of the current project. In this paper, we 
present our first results from a shallow survey over a large vol¬ 
ume, specifically targeting the brightest Ly a emitters. We are 
thus not expecting to find large numbers of candidates, but rather 
we are interested in their H a to Ly a properties and use ancillary 
information in well studied fields to connect these LAE candi¬ 
dates to other galaxy populations. 

We target well studied extragalactic fields to maximize the 
available multi-wavelength data, including observations with the 
Hubble Space Telescope (HST) in the ultraviolet, in order to 
add the third prong to our trident. In addition to the GOODS- 
North observations presented here, we are currently observing 
with NOT/MOSCA two Frontier Fields and Abell 1689 - lens- 
ing clusters with existing deep HST UV observations - using 
several new optical filters with deeper exposures than those pre¬ 
sented here (~ 10 hours per narrowband pointing). The total vol¬ 
umes probed per cluster will be similar to the central pointing 
of GOODS-North presented here, but reaching fainter emitters, 
also possibly magnified by the foreground clusters. We have 
started gathering Ha observations of Abell 1689 in the near- 
infrared, and we plan to observe our two Frontier Fields with 
CFHT/WIRCam in the near future. Once these observations are 
completed, they will be the subject of future Trident project pa¬ 
pers. 

Observing fields with HST UV coverage allows us to inves¬ 
tigate the leakage of Ly C photons from the galaxies we de¬ 
tect in the survey. These ionizing photons are intimately linked 
to the Ly a emission, as they are required for Ly a to form. 
The high absorption cross-section in H i for both types of emis¬ 
sion may cause them to escape together, probably through ion¬ 
ized cones or channels with low column densities and/or cov¬ 
ering fract ions, as mode l s of p orous mu lti-pha s e media sug¬ 
gest (e.g. IClarke & OevI l2(m iNeufeldl IT^ |L aursen et al. 
120131: IVerhamme et al.ll2014h . Observations and models of in¬ 
termediate redshift galaxies both support this picture, as LAE 
galaxies appear to emit m ore Ly C photons than what Lyman 
break galaxies (LBG s) do (llwata et al.ll2009HNestor et al.ll2013l: 
iMostardi et al.ll2M3h . Models also suggest t hat escape fractions 
of Ly C are hig her in lower-mass haloes (lYaiima et al.l 1201 H 
iMitra et ^l2013h . Note that all LAEs must show a Lyman break, 
and that selection effects will make an LAE survey find relatively 
more low-mass galaxies on average for a given exposure time. 

Adding t o this the observ ed steep UV LFs of faint galax¬ 
ies at z > 2 (lAlavi et al.ll2014t) . it therefore seems possible that 
low-mass Ly a emitting galaxies did have larg e enough Ly C 
escape fract i ons to ionize the univ erse (see e.g. iBouwens et aU 
I 2 OO 9 I l2012t iRobertson et al.ll2013h . and recent observations of 
the faint end of th e Ly a and UV LFs a t z ~ 6 - 8 are wel l in line 
with this picture (iDressler et al.l 120141: lAtek et al.ll20l'5l) . How¬ 
ever, due to the rapidly increasing opacity of the IGM beyond 
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Fig. 1. Fiducial map of pointings for MOSCA with filter NB401_6. 
Numbers correspond to pointings as described in Table [2] Dashed lines 
outline the GOODS ACS (red) and the CANDELS-Deep (blue) foot¬ 
prints. 

z > 4 dlnoue & Iwatall2()()^ llnoue et ^l2014l) . this Ly C escape 
can never be directly observed. Linking the H a, Ly a and Ly C 
emission together at z ~ 2 is therefore one of very few ways of 
quantitatively studying Ly C escape in galaxies. 

3. Observations and data reduction 

3.1. Custom filters 

The general idea of hnding a galaxy in a narrowband emission 
line survey is to do photometry of a large piece of the sky, us¬ 
ing a narrowband biter where the wavelength is centered on a 
redshifted emission line. When comparing the brightness of an 
object in this narrowband biter to a nearby broadband biter, any¬ 
thing with a strong emission feature in the narrowband biter 
will show a narrowband color excess (see e.g. lGuaita et al.ll2010t 
iNakaiima et al.ll2()T^ ISobral et al.ll2013L for recent descriptions 
of the narrowband technique and how it can be used). This there¬ 
fore gives a signature that can be used to bnd candidates for 
follow-up studies. 

One can of course combine narrowband biters at different 
wavelengths in order to see different emission lines. In our case, 
we want to look for galaxies not only in Ly a but also in H a, 
in a way that the redshift range is the same for both surveys. We 
therefore attempted to match near-lR narrowband biters with op¬ 
tical biters in terms of the redshifts of these emission lines. To 
this end, we had two custom biters made to exactly match ex¬ 
isting near-lR biters. We have also used three other, pre-existing 
optical biters, which extends our Ly a redshift range and thus in¬ 
creases the volume probed. We refer to the narrowband biters as 
e.g. NB392_6, where 392 and 6 are the central wavelength and 
width, in nanometers. NB392_6 and NB401_6 were designed to 
match H 2 and Brackett y near-IR biters. The properties of all 
the narrowband biters, and how they give an almost complete 
redshift coverage between z = 2.2 and 2.46, are summarized in 
Table □ 

Our Ly a observations were performed at the 2.56 me¬ 
ter Nordic Optical Telescope (NOT) on La Palma, Spain, us¬ 
ing the MOSCA instrument. MOSCA is a mosaic camera that 
is mounted at the Cassegrain focus, and it consists of four 
CCDs laid out to form a roughly 7.5 arcmin wide held of view. 


Table 1. Narrowband biter properties 


Eilter 


ID 

4c (A) 

4„ (A) 

z range 

ZP 

NB392 

6 

#134 

3917 

63 

2.20 - 2.25 

22.70 

NB401_ 

6 

#133 

4013 

64 

2.27 - 2.33 

22.72 

NB406_ 

6 

#121 

4058 

58 

2.31 -2.36 

22.52 

NB410 

6 

#131 

4102 

60 

2.35 - 2.40 

22.37 

NB412 

18 

#115 

4120 

180 

2.32 - 2.46 

23.66 


Notes. Nordic Optical Telescope ID number, central wavelength, width 
and corresponding approximate redshift range of Ly a for the narrow- 
band biters used in this study. The bnal column lists the average ze- 
ropoint for each biter, in AB magnitudes, acquired with a ~ 6 arcsec 
aperture. The uncertainty in all zeropoints is roughly 0.05 mag. 



Fig. 2. Transmission curves for the narrowband NOT biters used in this 
work (left panel) and how they relate to broad U, B and g’ biters (right 
panel). Note that the HST F435W biter is scaled to the entire instrument 
throughput. 

MOSCA’s CCDs are especially sensitive in the blue (with a 
quantum efficiency near 97% at 430 nm), and with few internal 
surfaces it is an ideal insbument for our purposes. 

We observed GOODS-North, using a central pointing of 
12h36”54% h- 62°14'50" (J2000). For NB401_6, we observed a 
grid of positions shown in Fig. [T] covering a surface area of ~ 

400 square arcminutes. Table |2]lists the exposure times and see¬ 
ing conditions for each pointing. The other narrowband biters 
where centered on the same coordinates but rotated 45 degrees to 
better match with the HST/ACS data footprint. The correspond¬ 
ing volume that we cover in the NB401_6 biter is roughly 83 
000 comoving Mpc^, with an additional ~ 60000 Mpc^ in sup¬ 
pleme nting biters. This is comparable to e.g. the iGuaita et al.l 
EolB MUSYC study, which covers ~ 135 000 Mpc^ of LAEs 
at z ~ 2.07, but goes considerably deeper. 

The Ha data for GOODS-North cover circa 22 arcmin^ 
(roughly the size of the view in Figure [TJ in a narrowband bi¬ 
ter matching the NB401_6 biter in redshift space. The details of 
these observations are described in a forthcoming paper. 

3.2. Data reduction 

The NOT/MOSCA data discussed in this paper were taken in 
March and May 2011 and February and May 2012. Table|2]sum- 
marizes these observations. 

Many of the basic reduction steps were developed from 
Tapio Pursimo’s notefl on reducing MOSCA data using the 
IRAIQmSCRED package. 

* http://WWW.not.iac.es/instruments/mosaic/reductions/reducti' 
^ IRAF is distributed by the National Optical Astronomy Ohservato- 
ries, which are operated by the Association of Universities for Research 
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Table 2. Observing log 


Date 

Filter 

Pointing 

Exp. time 

Seeing 

2011-03-01 

NB406_6 

1 * 

2.5h 

079 - 171 

2011-03-01 

NB410_6 

1 * 

1.5h 

171 - 172 

2011-03-02 

NB412 18 

1 * 

1.5h 

078 - 170 

2011-03-02 

NB406_6 

1 * 

1.5h 

078 - 079 

2011-03-02 

NB410_6 

1 * 

1.5h 

077 - 078 

2011-05-02 

NB401 6 

1 

2.5h 

078 - 079 

2011-05-03 

NB392_6 

1 * 

3.5h 

077 - 078 

2011-05-04 

NB401_6 

1,2,3 

2h,2h,0.5h 

078 - 170 

2011-05-05 

NB401 6 

2,3 

0.5h,2h 

078 - 171 

2012-01-27 

NB401_6 

5 

3h 

079 - 171 

2012-02-23 

NB401_6 

4 

4h 

079 - 170 

2012-05-20 

NB401 6 

2 

2 h 

173 - 175 

2012-05-21 

NB401_6 

3,5 

1.5h,lh 

079 - 175 

2012-05-22 

NB401 6 

6,7 

2 h,2h 

077 - 078 


Notes. Log of the NOT/MOSCA Observations. The filters listed are 
only the narrowband filters used in observations; complementing broad 
band observations (usually in the SDSS g’ band) where always per¬ 
formed. The "pointing" column corresponds to the numbers in Fig. [T] 
An asterisk is used to mark if the pointing was observed with a 45 de¬ 
gree field rotation. 


After the basic bias, flat and dark reduction steps in MS- 
CRED, the IRAF FIXPIX task was used to apply a bad pixel 
mask to the data. Next, the MSCSETWCS task from MSCRED 
was used to correct for individual rotation and distorsions in the 
four CCDs of MOSCA by changing the WCS projection. This 
step uses a plate solution which we constructed with the MSCT- 
PEAK task. 

The sky level was th en subtracted using Source Extractor 
(iBertin & Arnoutslll996h . and the images were combined to a 
single fits extension using the MSCIMAGE task. Cosmic ray 
cleaning was then applied using vO.4 of cosmics.py, a python 
interface to L.A.CosmicQ. The MSCIMATCH task was used to 
match the intensity between the frames, taking into account dif¬ 
ferences in airmass and seeing conditions, etc. The WCS of each 
frame was matched against the USNO-A2.0 star catalogue with 
the MSCCMATCH task, which was set to shift the WCS to the 
correct position. The science data were then stacked using the 
MSCSTACK task, which uses the WCS of the frames to align 
them before combining. 

The final astrometric solution was obtained by run¬ 
ning MSCCMATCH with a set of relatively bright HST b- 
band sources, in orde r to match our astrometry to that of 
iGiavalisco et alJ (l2004l) which is commonly used for sources in 
GOODS-North. We then used GREGISTER to register all of our 
ground-based data in all filters to a common set of physical im¬ 
age coordinates. This is a crucial step for doing the photometry 
with SExtractor in dual mode, where one image is used for de¬ 
tection and the other for photometry. 

For our ground-based data, we use aperture corrected pho¬ 
tometry where the correction is calculated from isolated bright 
sources in the field. We use aperture diameter sizes of 22 and 
60 pixels, corresponding to radii of 2.3 and 6.3 arcsec. De¬ 
tections were made using DETECT_THRESH =1.2 and DE- 
TECT_MINAREA = 32 in SExtractor, meaning that anything 
with 32 contiguous pixels that are all 1.2 times brighter than the 


in Astronomy, Inc., under cooperative agreement with the National Sci¬ 
ence Foundation. 

^ http: //obswww. unige. ch/~tewes/cosmics_dot_py/ 


background is considered a detection. This can be compared to 
-the seeing disk which for a full width at half maximum of 1 arc- 
sec would correspond to an area of about 17 pixels. Defining the 
thresholds this way allows for low surface brightness detections 
with unknown morphologies. 

We estimate the depth of our NB401_6 imaging to reach 24.3 
(24.6) AB mag at the 90% (50%) completion level, for the ar¬ 
eas with a total exposure time of ~ 4.5 hours. This corresponds 
to a limiting Ly a line flux (assuming a 20 A EW cut in the 
rest-frame) of 9.8 (7.3) x lO^'^erg s“*cm“^. In the NE and SW 
outskirts of the NB401_6 image, and in our other filters the ef¬ 
fective exposure time is ~ 2-3 hours, and the depth reached con¬ 
sequently 0.2-0.4 magnitudes shallower. 

3.3. Standard stars 

-Our custom made filters require the use of spectrophotometric 
standard stars f or the photo metric calibration. We used white 
dwarf stars from lOkj (Il990l) . as these stars do not show a strong 
Balmer break around 4000 A. We corrected the ab¬ 

solute fluxes according to lColina & Bohlin|(ll994l usually a shift 
on the order of -1-0.05 AB magnitudes). We then convolved the 
standard spectra with our measured narrowband transmission 
curves and calculated the AB magnitude of each standard star 
in each of our filters. 

The CFHT/WIRCam/WC8305 narrowband image was re¬ 
duced using the Simple-Wircam pipeline developed by Wei-Hao 
WangB NB401_ _6 is designed to match this filter, which was 
used for the matching H a survey, in redshift space. The H a re¬ 
sults will be discussed in a forthcoming paper, but the postage 
stamps in Fig |3] show a cutout for reference purposes here. 

3.4. Archival and catalog data 

A large set of multi-wavelength data already exists for various 
regions around GOODS-North. We have chosen to make use of 
ground-based UBVRIz broadband imaging using th e KPNO 4.0 
meter telescope and Subaru 8.3 meter, described in lCanak et al.l 
(l2004l) . We use the R-band detection catalog, using auto photom¬ 
etry. We chose this photometric catalog over the HST/ACS cata¬ 
log for several reasons. The covered area on the sky is larger, and 
includes the entire region that we examined with the NB401_6 
filter. The resolution is very similar to our data, as they are all 
ground-based, which makes the photometry easier to handle in a 
consistent manner. The KPNO U-band data is also important for 
our narrowband selection criteria. 

For our near-i nfrared data , we use the deep and 
Spitzer catalog from lWang et al.l (l2010h . who carefully matched 
Spitzer/IRAC data to each Ks source. To fill in the J and H 
bands where availab le, we make use of th e MOIRCS Deep Sur¬ 
vey near-IR catalog (iKaiisawa et alJ20T(ih . and calculate a J-Ks 
and H - Ks color. The MOIRCS Survey covers a horizontal strip 
through a large portion of the field. 

We do see a shift between the two Ks catalogs of about 
0.2±0.1 magnitudes. We investigated whether one of the two sur¬ 
veys would give a better optical-to-infrared color by measuring 
the average and median colors of all sources, but the end result is 
somewhere in between the two with quite a large uncertainty. We 
have therefore chosen to use the average value of the catalogs, 
and have added in quadrature an uncertainty of 0.2 magnitudes to 
all of our near-infrared and Spitzer data points. This reflects our 
uncertainty in the relative optical to near-infrared photometry. 

http: //WWW. asiaa. sinica. edu. tw/~whwang/idl/SIMPLE/ 
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Fig. 3. Postage stamps for an example candidate galaxy. The field of view is 10 x 10 arcsec. The large RGB picture shows an HST bvi color 
composite made with the scale in ea ch color channel rel ated to the true total intensity, with an arcsinh intensity scaling relation. This coloring 
scheme is similar to that described in lLupton et alj j2004l) . The white circle has a one arcsec radius - typically twice the size of the ground-based 
seeing disk. Postage stamps for all candidates are compiled in the Appendix. 
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Fig. 4. Narrowband selection for our NB401_6 filter. Dotted lines show 
the cut-off. Photometry for all data is indicated by small gray points. 
Large blue points indicate candidates that meet both the excess and the 
BzK criteria. The blue data point on the top left shows the typical size 
of the uncertainty. 



Fig. 5. BzK selection for candidates in the NB401_6 filter. T he dotted 
line shows the cut-off, corresponding to the lDaddi et alJ(l2004h BzK cri¬ 
terion. Photometry for all data is indicated by small gray points. Large 
blue points indicate candidates that meet both the narrowband color se¬ 
lection and BzK criteria. The blue data point on the bottom right shows 
the typical size of the uncertainty. 


4. Candidate selection 


There are several possible approaches to finding LAE candidates 
in a narrowband survey. A common way is to simply measure 
a narrowband to broadband color, where the broadband covers 
the nearby continuum. If a strong emission line falls within the 
narrow filter, this will produce a narrowband excess and possible 
candidates will therefore show a color signature. This method 
is relatively stable if the broadband is centered approximately 
on the narrowband, or if the continuum is expected to be flat 
throughout the probed wavelength range. 


In this simplified case, the color signature in the narrowband 
filter can easily be estimated. The flux density in the narrowband 
filter can be described as 


fv,NB - fyNB + /v 


fcont 

fv,NB 


f fr'(X}f2TNB(X}e-^'^^(^'^dA 
f ■^TisiB(X)dA 


( 1 ) 


Where and are the flux density contributions from 
the emission line and continuum, Tnb is the transmission curve 
of the filter, and the IGM transmissivity is described as a sim¬ 
ple exponential with appropriate t/cm fe.g. lMad^ll995h . Such 
IGM absorption bluewards of Ly a can be significant, especially 
at higher redshifts. The expression for the broadband filter will 
have an identical appearance. 

The continuum can be described as a power law with slope 
a, so that f™"' = —A°‘, where A is a normalization factor. If we 
assume the Ly a line to be a delta function, and the IGM effects 
to be negligible, we can simplify (HJ to 
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FelTnb fJVA--^Tf,BU)dA 
f ^TjvB(A)dA f -^Tji/b (A)dA 


where Fee is the flux in the emission line and T/^b is the 
integral averaged transmissivity of the Alter. Describing Fee in 
terms of the equivalent width, Fee - EW x NA°‘~'^, we arrive 
(with (|2|) at an expression for the flux in a given Alter as a func¬ 
tion of equivalent width and the continuum parameters N and a. 
The color for a given narrowband to broadband is then simply 
given by the usual magnitude definition. 


niNB - niBB = -2.5 logjo 



(3) 


With a flat spectrum in /y (a - 0) and an observed EW cut 
of ~ 65 A (20 A in the rest-frame at z ~ 2.3), (O and the method 
described above would give us an excess with a color cut g' - 
NB ^ 0.6, or t/ - AB 0.7. 

Note that the NB412_18 Alter is three times wider than the 
other filters. This means it is probing a three times wider red- 
shift space but at the same time reaches a lower significance due 
to decreased contrast between continuum and line emission. We 
can only reliably resolve sources in this Alter with an equiva¬ 
lent width in the observed frame similar to the width of the Alter 
(>180 A). 

For a more realistic calculation, we need to take the com¬ 
plications mentioned above into account. Our narrowbands have 
wavelengths near the edges of our broadband Alters, and more¬ 
over, the continuum on either side of Ly a can not be expected 
to be flat. Absorption of Ly a in H i gas in the IGM will tend to 
attenuate the spectrum bluewards of Ly a. A strong continuum 
slope can therefore create false candidates, or make real sources 
disa ppear, if a classic narr owband excess technique is used (see 
e.g. IHaves & Ostli'nll2006l) . We have therefore chosen to com¬ 
bine two broadbands, the U and g' bands, on either side of Ly a 
atz ~ 2. 

In order to estimate our color selection criteria, we first cre¬ 
ated model spectra using our SED fitting code GalMC (see Sec¬ 
tion ISll. We convolved these spectra with our measured Alter 
transmission curves to calculate the expected brightness in each 
filter. 

Using an EW(Ly a) cut of ~ 20 A (rest-frame), we start by 
a color cut U - NB > 0.7, corresponding to the limit of no IGM 
absorption. This color also has a significance criterion, so that 
the photometric error is not larger than the color itself. To al¬ 
low for reddening up to E{B - V) = 0.5, we set our limits at 
g' - NB > Q.\ and g' - NB > -1.2{U - NB) + 1.4. Allowing 
for stronger reddening could introduce a group of interlopers at 
slightly different redshift, that may show a strong U - NB ex¬ 
cess due to the Lyman break, but have no actual Ly a signal in 
the narrowband. Our limit is therefore chosen to balance these 
effects as well as possible. Many of our candidates have a large 
photometric uncertainty due to this being a shallow survey, and 
we therefore also chose this double-excess method to get a more 
reliable selection, although we will always have some interlop¬ 
ers in our candidate sample. The -1.2 slope was calculated from 
our GalMC models, and is set by the relative widths of the U and 
g' Alters. 

In order to deal with possible low redshift interlopers with 
[On] hlTl A emission, we combine these criteria with the stan¬ 
dard BzK = [z - K) - (B - z) > -0.2, designed to select star 


forming galaxies at z > 1.4. Figures]?] and |5] show these criteria 
applied to our NB401_6 data. Using the BzK criterion requires a 
detection in these broadbands, and this criterion reduces our to¬ 
tal number of LAE candidates from 42 to 30 sources. Given that 
the broadband data is considerably deeper than our narrowband 
data, and that we are looking specifically for bright emitters, we 
use this also as a criterion to remove spurious sources. 

We correlate our sa mple of candidates to the spectroscopic 
catalog ues compiled by iBarger et al.l (l2008h and lAdams et al.l 
d201 ih . These catalogs only contain sources for which the red¬ 
shift could be determined with several spectral features. We And 
flve candidates to already have spectra consistent with the ex¬ 
pected narrowband wavelengths (as indicated in Tables 01 - |6]l, 
but there are also spectroscopically discovered LAEs which are 
too faint to be detected in our sample. A thorough analysis of the 
Ly a and H a brightness of galaxies which are already spectro¬ 
scopically conflrmed will be included in a forthcoming paper. 

We match o ur candidates to the Ch andra Deep Field North 
X-ray survey ( Alexander et alJ l2003h and deep VLA data 
(iMoirison et alJ 201 Oh in order to remove AGN. Anything with 
an X-ray luminosity > 10"^^ erg/s or with a 20 cm detection 
> 10^"^ W/Hz is considered an AGN. We remove 5 candidates 
based on these criteria, leaving 25 LAE candidates for further 
analysis. Our AGN fraction of ~ 17% is broadly consistent with 
similar works (see Section lSTI) . 


5. Warm dust properties 


We have performed a match again s t the p hoto metric GOODS- 
Hersch el catalogs of lElbaz et ^ d201 ih and iMagnelli et al.l 
( 2013h. Th e latter catalog combines the deeper data from 
Elbaz et ^ (1201 Ih with the PEP 100 and 160 pm data of 


Lutz et al. ( 201 ih . 19 of our candidates are in the held of view of 


these surveys. 

We And flve candidates with 24 pm Spitzer/MIPS detec¬ 
tions within 1 arcsecond of the catalog coordinates. These are 
NB401_6-10,-14,-19 and NB392 6-l,-2. Three of thes e flve 
candidates have detections in the IMagnelli et al.l (1201 3h PEP- 
Herschel catalog; NB401_6-19, NB392_6-1 and NB392_6-2 
(see SectionIQ. 

Interestingly , NB4 01 6-10 has a detection at 160 pm in 
the lElbaz et SJ (1201 Ih catalog, but with no detection at 100 
/rm and not be i ng inc luded in the 3cr catalog of the combined 
IMagnelli et al.l (1201 3h data. It seems that the original flux may 
have been overestimated, and that this candidate no longer meets 
the significance cut in the com bined catalog. At th e same time, 
NB392_6-1 is not listed in the lElbaz et al. (1201 ih catalog, but 
with the deeper data it is det ected at 3cr. _ 

For the three matches in IMagnelli et al.l (1201 3h . we have fit¬ 
ted an M82 template to the far-IR data points, i.e. only to the 
datapoints from 24 pm or longer wavelengths, to derive an es¬ 
timate of the total IR luminosity, Ljb, integrated from the rest- 
frame 8-1000 pm. Figure |6| shows the fit to the data, with esti¬ 
mates of log(L/R/Lo) in the bottom r ight corner of each graph. 
Also shown are two AGN templates dPolletta et al.ll200^ - one 
quasar template and a Seyfert 2 galaxy template - as well as an 
Arp220 template for comparison. 

Assuming the M82 template fit value as our value for Ljb, we 
compute the total IR to UV luminosity, IRX = log(L/R/Li6oo), 
where Ligoo is the monochromatic luminosity at rest-frame 1600 
A, cal culated from the photometry in the V band. iMeurer et al.l 
(II 999 II found a tight correlation between the rest-frame UV slope 
yS and IRX from a local sample of starburst galaxies, indicating 
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Fig. 6. Fits to the far-infrared data points from the lElbaz et alj (1201 Ih GOODS-Herschel catalog, using two templates of star forming galaxies (M 
82 and Arp 220 in solid red and dashed blue, respectively) as well as two AGN templates (QS02 and Seyfert 2 in dot-dashed green and dotted 
cyan). Note that only the infrared > 8yum catalog data were used in the fit, and the optical to near-infrared data are only shown here for comparison. 
The numbers in the bottom right of each panel show the integrated 8-1000 /rm luminosities in solar units for each template. 



Fig. 7. Left panel: p vs. IRX plot for our three cand idates detected in G OODS-Herschel. Note that the iMeurer et al J lll999h data and relation 
have been corr e cted w ith t he 1.75 fac t or as q uoted by [Reddy et alj ( 1201 Oh . to convert from Lfm to the total Lm. Also shown are the more recent 
iTakeuchi et alJ (12012h and iBuat et al.l (12012h j8-IRX relations, derived for local galaxies and z ~ 1 - 2 so urces, respectively . Right panel: IRX 
vs. bolometric luminosity for the same candid ates, compared to th e sample of BzK selected galaxies from iDaddi et al.l ( 120071) and the empirical 
relation for BM/BX-selected galaxies found bv iReddv et al.l(l2006h . 


that the energy absorbed by dust in the UV is being proportion - 
ately re-emitted at far infrared wavelengths. iReddv et al.l (l2010h 
found that this relation breaks for ultra-luminous infrared galax¬ 
ies (ULIRGs), probably due to saturation as part of the UV-bright 
star forming regions are being completely obscured by dust. 


We estimated the fi slopes for our three candidates by fitting 
a power-law through the observed V, R and I magnitudes, cor¬ 
responding roughly to the 1600-2600 A range in the rest frame. 
Overplotting our three candidates in the left panel of Figure |7] 
they appear to fall into the sam e locus as the z ~ 2 ULIRG 
sample from iReddv et al.l (1201 Oh . Indeed, we find Ljr > lO'^L© 
for two of our three candidates, which is t he very definition of 
a ULIRG. In this panel we also show the iMeurer et al.l (Il999h 
relation as corrected by iReddv et al.l (l2OI0l) for comparison . 


Note that the GALEX IRX ratios used byjMeurer_etalJ( 1999) 


are believed to have been overestimated (lOverzier et al.l 


2011 


ITakeuchi et al.ll20l^ . We show for comp arison the new r elatio n 
derived for local galaxies as presented bv lTakeuchi et al.l (l2012l) . 


and th e relation derived for z ~ 1 - 2 sources from iBuat et akl 
(I 2 OI 2 I) . In the right panel of Figure|7]we see how our candidates 
com pare with the IRX- lumi nosity ratios of z ~ 2-3 sources 
from iDaddi et al.l (l2007l) and IReddv et al.l (l2006l) . All three can¬ 
didates appear to lie within the normal spread, indicating that our 
FIR-detected LAEs do not possess unusual UV to IR properties, 
compared with other samples at similar redshift. 

We estimate the energy budget for our candidate LAEs by 
calculating the amount of flux absorbed by dust in the optical- 
UV part of the spectrum, and express it in terms of luminosity 
as Lahs,uv- We create model spectra based on simple SED fit¬ 
ting (modelling only age and E(B-V) for constant star formation 
rate; see model A in Section |6]l, both with the best fit value of 
E(B-V) and with a value of 0. We then integrate both spectra 
for each candidate, and estimate the amount of flux absorbed by 
dust as the difference between the two spectra. Note that Labs,uv 
is not an observed quantity but an estimate of the amount of dust 
absorption, expressed in terms of the absorbed luminosity. The 
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Fig. 8. Far-IR luminosity versus dust-absorbed luminosity as estimated 
from the E(B-V) values derived from SED model A. Our three far-IR 
detections are marked with magenta stars. Blue triangles show estimates 
at our approximate detection threshold in the far-IR of Lm ~ 10**'*Lo. 
Black triangles show upper limits to the dust absorbed luminosities, and 
have been shifted by -0.02 dex along the y-axis for clarity. The dotted 
black line shows the 1:1 relation. 
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results are shown in Fig. [8j plotted against the estimated far-IR 
luminosity. Our three (U)LIRGs are found among the candidates 
that appear to have the highest amount of Lats.uv, as expected. 
All candidates show Labs.uv < Ljk, which is also expected in 
the scenario where a dust screen model may underestimate the 
total dust content. Since these estimates of Lahs,uv should cor¬ 
respond to a lower limit to L/s, we thus expect at least one of 
them to be bright enough to be detected as a ULIRG in the far- 
IR (NB401_6-19, which is also our brightest candidate in the 
far-IR). A handful of the rest of the candidates would appear to 
be at least in the LIRG regime (11 < \og(Li;v,abs) < 12), and as 
discussed, two of these are also detected as LIRGs or ULIRGs 
in the far-IR. 


5. 1. Comparison to previous work 

Our findings are broadly consistent with several previous works. 
Although the LAE population is normally associated with small, 
dust-free systems, evidence of Ly a emission from large, dusty 
galaxies does exist in the literature, but the AGN contribution to 
the Ly a emissivity is not always well known. 

IChanman et al.l (l2005h present a spectroscopic follow-up 
study of 98 sub-millimeter galaxies detected at 850 pm and deep 
VLA 1.4 GHz data, and obtained reliable redshift estimates for 
73 of these. 25 of those sources show a "strong" Ly a line. The 
AGN fraction is not certain, as they cannot rule out faint AGN 
features, but they note that the Ly a to [C iv] A 1549 line ratios 
are generally consistent with what you would expect from nor¬ 
mal starbursts. They conclude that a patchy or inhomogenous 
dust distribution might allow Ly a to escape from these systems. 

ICasev et"^ (1201 2h performed a spectroscopic study of 767 
galaxies individually detected in Herschel/SPIRE data. 36 of 
those galaxies are found to be at z > 2, with L/r ~2-6x1 0 '^Lq, 
the majority of which show Ly a in emission. 17 of the Ly a 
emitting galaxies do not show any [C iv] emission, and they es¬ 
timate an AGN fraction of roughly ~ 25 %. 

HST/COS spectroscopy in the ul traviolet of 11 _ local 

ULIRGs is presented in recent work by iMartin et al.l (l2015h . 8 
of the ULIRGs show Ly a in emission, three of which appear to 
be AGN dominated galaxies. Eor the five remanining ULIRGs 


they estimate AGN to contribute no more than 15% to 35% to 
the bolometric luminosity. 

lOteo et akl (1201 2al) perform a match of 56 spectroscopi¬ 
cally detected LAEs at z ~ 2-3.5 in GOODS-South to the 
GOODS/Herschel catalogs and find 10 of them to be detected 
by MIPS at 24 pm, four of which have detections in the PACS- 
160//m band within 2 arcseconds of the optical detection. All 
four candidates show integrated luminosities L/r > IO'^Lqj plac¬ 
ing them in the ULIRG regime, while the other six galaxies de¬ 
tected in MIPS are at least in the LIRG regime (L/r > IO'^Lq)- 
Their study of 23 LAEs at z ~ 0.3 with mid-IR/EIR detections 
showed only one such IR-lu minous galaxy, with the r est of the 
sample having Ln ? < 10**Ln jOteoet al.ll201 iL [2012bh . To echo 
the conclusions of Oteo et al.l(l2012ah . this indicates that the Ly a 
selection technique might trace different galaxies at different red- 
shifts; as seen for example in GALEX LAEs which are generally 
not as luminous a s their higher redshift counterparts (see e.g. 
ICowie et al.ll201 1|) . 


On the other hand, our results d o not appear to agree w ith 
a similar study recently performed bv iKusakabe et al.l (1201 5h on 
LAEs at z=2.18 in the GOODS-South field. After stacking 218 
LAEs in the Herschel bands they report only an upper 3cr limit 
on L/r of 1.1 X IO'^Lq. Our study does not include the same 
removal of 24 pm detections and may therefore be more vulner¬ 
able to confusion from nearby sources. However, judging from 
the multi-wavelength imaging we have available, our sources are 
well-isolated (with the possible exception of NB401_6-19) in the 
Spitzer/IRAC imaging, which was used to identify detectio ns for 
the Herschel sample. According to lKusakabe et akl (1201 5h . none 
of the 218 LAEs in the GOODS-South sample were detected 
in MIPS 2Apm imaging at the 3cr level. As previously stated, 
out of our 19 candidates overlapping with the MIPS footprint, 
five have a counterpart in t he 24p m band within 1 arcsec, as re- 
ported by MagjellLet^l (120131) . The 218 LAEs discussed by 
iKusakabe et al.l ( 1201 5l) are however considerably fainter on av¬ 
erage than the sample presented here (Table[3]i; the average Ly a 
luminosity is 6.6x 10"^' erg/s, and the maximum is only 6.6x 10^^ 
erg/s (Kusakabe 2014, private communication). The average lu¬ 
minosity for our sample is 4.2 x 10"^^ erg/s (see TableO. 


IWardlow et al.l (1201 4h stacked three samples of 126, 280 and 
92 LAEs in the Extended Chandra Deep Eield-South (ECDE- 
S) at z~2.8, 3.1 and 4.5 respectively in Herschel/SPIRE 250, 
350 and 500 jUm i mages as well as 870 pm data from LABOCA 
(IWeiB et al.ll2009l) . Similarly, they put Icr upper limits on L/r of 
5 - 8 X IO'^’L©, corresponding to roughly 2 - 3 x lO^L© at 3cr 
for all three redshift ranges. Only one of our candidates is bright 
enough to be detected at these longer wavelengths; NB401_6- 
19 (at 250 pm), which is also the brightest in the far-IR of our 
candidates. 


In light of these somewhat conflicting results, it is clear that 
large samples of bright LAEs, in particular at the strong star¬ 
forming epoch around z ~ 2, with matching far-IR surveys could 
test if this is a cosmic variance issue. 


As a final note, a preliminary analysis of our 
CEHT/WIRCam data of Ha candidates matching the red¬ 
shift of our NB401_6 Ly a data finds both NB401_6-19 and 
NB401_6-10 also show significant Ha emission, which would 
be expected for bright infrared galaxies with high star formation 
rates, if the star forming regions are not too obscured by dust. 
The details of the full H a to Ly a match will be discussed in a 
forthcoming paper. 
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6. SED fitting 

We co mbin e catalog data f r omlCapak et al.l (l2004l) . IWang et ^ 
(1201 oh and iKaiisawa et alJ (1201 Oh as described in Section 13.41 
In order to perform SED fitting on a relatively homogeneous 
and self-consistent data set, we chose to use archival Sub- 
aru-i-KPNO UBVRIz data rather than HST/ACS data. A collec¬ 
tion of "postage stamps" from our imaging data is presented as 
an example in Figured The full set is included in the Appendix. 

To perfor m the SED fitting, w e use version 2.0 of the 
GalMC code (lAcauaviva et al.ll201 Ih . with GALAXEV synthe¬ 
sis modeE_of_S^_C^rlo^& G. Bruzual (2007, updated version 
of iBruzual & Charl^l200.3h ') and a Salpeter initial mass function 
(IMF). We also ran some tests, comparing the original (2003) 
and the more recent (2007) Charlot-Bruzual model, as well as 
Chabrier and Kroupa IMFs, but found no strong differences for 
the galaxies dominated by yo ung stellar population s that we are 
mainly interested in. We use a lCalzetti et al.l (l2000h dust absorp¬ 
tion model. For an overview on these different models and the 
general effects they ha ve on SED fitting, see the discussion in 
lAcauaviva et alJ (1201 Ih . GalMC uses the information on the re¬ 
processing of ionizing photons in order to add nebular emissio n 
lines, similar to the approach of ISchaerer & de BarrosI (l2009h . 
The Monte Carlo Markov chains were run using 20 000 samples. 
The output chains were then analyzed using the April 2013 ver¬ 
sion o f the GetDist program from CosmoMC (iLewis & Bridlel 
[ 2001 . 

For all of our SED fitting we have chosen to remove the U- 
band, corresponding to the spectrum blueward of Ly a, from the 
fit. The co de does includ e a correction for IGM absorption ac¬ 
cording to lMadai](ll995h . but this is only an average correction. 
Since this parameter is not known for the individual galaxy, in¬ 
cluding the U-band could skew the results. We have also chosen 
to subtract the contribution of the Ly a flux to the B-band. This 
is done simply by assuming the contribution from the continuum 
to be the same in B as in the narrowband. The uncertainty in the 
B-band flux was then increased by adding in quadrature the pho¬ 
tometric uncertainty of the narrowband flux contribution. 

The metallicity is a notoriously difficult parameter to con¬ 
strain in SED fitting, as it has a relatively modest effect on the 
overall SED shape and is also strongly degenerate with the dust 
reddening. We therefore chose to use a fixed metallicity of 0.2 
Zq and only allowed the reddening to vary. We held the redshift 
fixed at the expected value, corresponding to Ly a being centered 
in the narrowband filter. 


6.1. Selection of SED fitting parameters 

Galaxies that we observe at z ~ 2 have already lived through over 
2 Gyr of cosmic time. In that time, these galaxies have likely 
gone through several merging events with epochs of heavy star 
formation as well as times of more quiescent evolution. What we 
observe at this moment represents a snapshot in time somewhere 
along that evolutionary path. To accurately resolve the entiry his¬ 
tory of such a complex system based solely on the spectral en¬ 
ergy distribution is a challenging task, which comes down to a 
matter of statistics and averages. 

Based on certain assumptions, we can see which parameters 
would fit the observed SED the best. Rather than only assuming 
the standard SED fitting parameters and simply compiling our 
results, we wanted to see if we could say anything more defini¬ 
tive about the star formation histories of our LAE candidates. We 
therefore decided to run SED fitting using three different setups: 


- Model A corresponds to the common approach of fitting for 
the age, reddening and mass of the galaxy, using a constant 
star formation rate. For this model, the age values represent a 
form of luminosity-weighted age, based on the entire optical 
to near-IR SED. This age is in other words difficult to inter¬ 
pret in a physically meaningful way. A high (low) age is usu¬ 
ally associated with a red (blue) spectrum, but this quantity 
is of course strongly degenerate with the dust reddening. In 
order to attempt to break some of these degeneracies, and ex¬ 
plore the parameter space in more detail, we also attempted 
two more unconvential models. 

- Model B assumes an age of the galaxy equal to 2.2 Gyr, cor¬ 
responding roughly to the age of the first stars in the uni¬ 
verse at z ~ 2. This is not a typical value found for the stellar 
model age of LAEs, but the idea of this model is to see if 
we can constrain a general trend in the star formation history 
throughout the formation of the galaxy. The star formation 
rate rj) is then allowed to change according to a simple expo¬ 
nential law; (j) - 1 /|t| exp(-f/T). By allowing the r parameter 
to take any value |t| > 1 Myr, the star formation may either 
decrease or increase with time, corresponding to positive or 
negative r values, respectively. The exponential models are 
restricted so that |age/T| < 13.78. 

- Model C assumes two separate populations, each with con¬ 
stant star formation histories, where one is fixed to an age 
of 2.2 Gyr and the other is restricted to be younger than 10 
Myr. This age corresponds roughly to the maximum age for 
which a simple stellar population can be expected to produce 
ionizing continuum photons giving rise to nebular lines, in¬ 
cluding Ly a. This model is thus similar to the idea in model 
B, where we want to see how much more recent star forma¬ 
tion is needed in order to reproduce the observed spectrum. 
We can then estimate the relative contributions of the two 
populations to the total mass as a measurement of whether 
the old or young stars dominate the SED. 

The idea of using these three models is not to try to find the 
one model that best represents the spectral shape of our data. In 
fact, all three models can usually find reasonable fits to all of 
the observed photometry. None of the models represent a truly 
realistic view of the star formation history and so are not meant 
to be compared in that sense. Instead, our idea is to try to see 
what SED modelling can tell us about the star formation rates in 
general, under certain assumptions. 

6.2. SED fitting results 

The results of our fits for all our LAE candidates are presented in 
Tables |4]|5] and |6]for models A, B and C respectively. The values 
given as best fit parameters are those corresponding to the low¬ 
est found. However, the probability density space may show 
a complicated structure, and the following "low" and "high" 
columns listed after each parameter represent the 68% confi¬ 
dence intervals marginalized over the entire parameter space. 
Note that the least-squares solution may even in some cases be 
located outside this region. This represents a case where most 
reasonably good models favour certain values for that parameter 
(e.g., a low age), while the single model with the lowest hap¬ 
pens to correspond to a different value for that parameter in that 
particular case (e.g. a high age). The marginalized limits there¬ 
fore give a slightly broader picture of which parameter values 

can still yield acceptable fits^ _ 

The R-1 value shows the lGelman & Rubinl(ll992ll R statistic, 
which is a measurement of how similar the chains are towards 
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Fig. 9. EWLya and LLya vs. Stellar mass, derived from SED fitting in model A. Downward pointing triangles denote upper Icr limits. 
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Fig. 10. Left: Age vs. Mass for SED model A, using constant star formation rate. Note that NB401_6-7 is not included in this graph as the age 
could not be constrained for this candidate. The mass in all three panels is the total mass that has gone into stars, i.e. the integral of the star 
formation rate over the age of the galaxy. Middle: Birth parameter vs. Mass for SED model B, using an exponential star formation rate (Eq.|4l(. 
Lower limits in the birth parameter correspond to lower limits in 1/t and are denoted by upward pointing triangles. Right: Birth parameter vs. 
Mass for SED model C. Note the logarithmic scale on the y-axis. Upper limits (Itr) are denoted with a downward pointing triangle. Four of the 25 
candidates are not shown in this panel, for which the mass of the old population could not be recovered (see text). The mass is the total mass of 
the two populations combined. 


the final stages of the fitting. Ideally, all chains will end up in the 
same region of parameter space, indicating a good convergence. 
In cases of strong degeneracies with multiple local minima, this 
number may be larger, indicating that the different chains favour 
different solutions. R-1 < 0.1 is typically considered as accept¬ 
able convergence. Any value above 0.1 is typeset in bold in Ta¬ 
bles |4l|5] and |6] for clarity. 

The general picture that emerges from the different models is 
not unexpected; forcing an old age generally shifts the extinction 
values down, making the spectrum less reddened due to dust. For 
the exponential star formation in model B, a negative r value is 
often favored, corresponding to exponentially rising star forma¬ 
tion, again filling in the blue spectrum by adding more newly 
formed stars. 

In model C, we wanted to test how important recent star for¬ 
mation truly is for reproducing the spectrum. Our results indi¬ 
cate that many of these LAEs seem heavily dominated by a re¬ 
cent burst of star formation, while some are consistent with a 
slow, steady rate. Fig. [TO] shows the age or the burst parameter b 


against the mass for our three models. The b parameter is defined 
as the current star formation rate (SFR) divided by the average 
past SFR. For model C, there are four candidates for which the 
best fit mass in the old population is negligible, corresponding to 
a system where a ~ 10 Myr population completely dominates the 
spectrum. These four systems are also clearly among the lowest 
mass systems in model C (see Table |6]l. Since we fixed the age 
of the old population to 2.2 Gyr, the resulting burst parameters 
are overly high (note the log scale in the right panel of Fig. fTOl i. 
and would for those four cases imply an essentially infinite b 
parameter. These four candidates are therefore not shown in the 
figure. Note that the b parameters derived for both model B and 
C are not to be interpreted literally, due to the heavily restricted 
mode of star formation in the two models. Rather, we show the 
b parameters here as a way to understand the ratio of young to 
old stars in our models, and see the results only as an indication 
that the majority of our LAE candidates show signs of recent 
elevated star formation. 
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For model C, the least massive systems are more inclined 
towards being dominated by a young population, and vice versa. 
When including all data points, the Pearson r coefficient in log 
space between the b parameter and the mass is -0.659, with a 
Spearman rank p - -0.576. There is a tendency for this also in 
model A (left panel of Fig. [TOl i. where the ages appear slightly 
lower on average for the less massive systems, although the r 
and p coefficients are only 0.261 and 0.303 for this case and 
the uncertainties are large. The trend is even less clear in model 
B (middle panel of Fig. fTOl). where the birthrate parameter is 
estimated as (e.g. lMadgwick et'ani2003h : 


b — 


SFR(f) 


(SFR(f)>past 


[t - tf) e-{t-t,)lT 
T 1 _ e-(!-tf)iT 


(4) 


In eq. m f - is the difference between the current time 
and the time when star formation first began (i.e., the age of the 
galaxy), and r is the characteristic time-scale of star formation. 
Note that although we can estimate the birthrate parameter us¬ 
ing eq. a an exponential SFR implies a smoothly increasing or 
decreasing rate, rather than a sudden burst. We can thus only use 
model B to say whether our candidates appear to have increased 
or decreased their star formation rate in a more general sense. 

The combined picture of these three models is clear; most, if 
not all, of our candidates show signatures consistent with recent, 
strong star formation, and the relative amount of young stars ap¬ 
pears to be higher in the least massive systems. This is broadly 
consistent with a scenario where our LAEs are galaxies that have 
recently undergone an episode of strong star formation, across a 
large spread of galaxy masses. 

The Ly a luminosities and equivalent widths are presented 
in Table[3] These values were calculated using the continuum fit 
around Ly a from model A. Comparing our sample to similar 
analysis of z ~ 2 LAEs in the literature, we can see for our 
model A, plotted in Eig|3 our sample c ontains relatively m assive 
(i.e. lu minous) galaxies. Co mpare e.g. lHagen et al. ( 2014), thei r 
Eig. 7, IVargas et al.l (|201^. their Eig. 5 or iGuaita et al.l ( 201 Ih . 
their Eig. 3. Nilsson et alJ (1201 ih also found only higher stellar 
masses, similar to ours, for a sample of LAEs compiled in a 35 
X 34 square arcminute survey using a two-meter class telescope. 

The average star formation rate derived from SED fitting of¬ 
ten shows signs of being elevated in samples of LAEs, as a func - 
tion of stellar mass (e.g. lVargas et aljF2014l : lHagen et al.ll2014l) . 
This is consistent with a picture in which small, starbursting 
galaxies produce large amounts of ionizing photons, which are 
then reprocessed as Ly a which escapes relatively easily from 
the lower optical depth of the smaller systems. Looking at the 
star formation rate against stellar mass in Eig.[TT] our candidates 
are equally consistent with BzK-selecte d samples of strongly 
star forming galaxies (iDaddi et al.l 2007|) as with L v a-only se¬ 
lected galaxies at z ~ 2 - 2.5 ( Whitaker et al.ll2()T^ . This is not 
surprising given our BzK selection criterion, as well as our aim 
to target brighter galaxies in this survey. 

Eigure [12] shows the relationship between the mass and dust 
reddening as found in our SED fitting model A. The data points 
are marked by a colored circle, representing the amount of flux 
in the Spitzer/MIPS 24 pm band (see Section|5]l, as indicated by 
the colorbar on the right hand side. Non-detections are white in 
the figure. All five candidates with 24 pm detections are located 
in the upper half (corresponding to the highest dust obscuration) 
of the 19 candidates overlapping with the MIPS footprint. This 
indicates a strong connection between the amount of reddening 
of the spectrum in the optical and ultraviolet to the emission at 


Table 3. Ly a luminosity and EW 


Candidate 

Tiv Q'(10'‘'^ergs^‘) 

EWiv a 

NB392_6-1 

4.2 + 3.3 

116 + 90 

NB392 6-2 

6.8+ 4.0 

190+112 


NB401_6-1 

4.2 + 4.2 

66 + 65 

NB401 6-2 

11.3 + 8.5 

160+120 

NB401_6-3 

2.3 + 2.4 

112+117 

NB401_6-4 

5.5 + 3.1 

47 + 26 

NB401 6-5 

2 .8+ 2.1 

93 + 71 

NB401_6-6 

3.1 +2.1 

297 + 203 

NB401_6-7 

2.8 + 3.7 

24 + 32 

NB401 6-8 

4.1 +3.0 

150+107 

NB401_6-9 

3.5 + 2.4 

75 + 52 

NB401 6-10 

1.9+ 2.1 

35 + 39 

NB401_6-11 

5.1 +4.3 

30 + 26 

NB401_6-12 

5.2+ 5.3 

87 + 89 

NB401 6-13 

1.5 + 1.5 

58 + 55 

NB401_6-14 

2.5 +2.1 

58 + 49 

NB401_6-15 

2.4 + 2.7 

66 + 74 

NB401 6-16 

5.6+ 5.0 

115 + 102 

NB401_6-17 

2.3 + 2.3 

98 + 99 

NB401_6-18 

1.9+ 2.3 

54 + 65 

NB401 6-19 

3.3 +2.2 

60 + 41 

NB401_6-20 

2.1 + 1.5 

125 + 87 

NB401 6-21 

7.6 + 6.5 

120+103 

NB406 6-1 

6.7 + 5.4 

142+114 

NB412 18-1 

7.6 ± 10.2 

103 + 139 


Notes. Ly a luminosities and equivalent widths for our LAE candidates. 
These quantities were calculated using the SED fit from model A for the 
continuum estimation. Candidates with detections in the far-IR (Sec- 
tion|5]( are marked in bold. 



Fig. 11. Stellar mass versus star formation rate for our LAEs, derived 
from the S ED fitting in model A. The relation found for BzK selected 
galaxies bv lDaddi et alJdTOOTI) (solid line) is also extrapolated to lower 
masses (dashed line). The IWhitaker et aP (1201 2f) relations for z = 2.0 
and 2.5 are shown in red and black dotted lines, respectively. 


infrared wavelengths. The 24 pm band corresponds to roughly 8 
pm in the rest-frame, which is where strong PAH emission and 
hot dust dominates the spectrum. 
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Fig. 12. Stellar mass vs. E(B- y), derived f rom SED fit, model A. The 24 
fim MIPS flux, as reported in lElbaz et al.l ( 1201 Ih . is indicated by color. 
Open circles represent no detection at 24 jim. 


7. Conclusions 


We have performed a survey of Ly a emitting galaxies (LAEs) 
at z ~ 2 in GOODS North using narrowband filters at 
NOT/MOSCA. Drawing upon the large quantity of ancillary data 
in the field, we use the entire optical to infrared catalogues to at¬ 
tempt to classify these LAE candidates. 

We perform a detailed SED analysis using GalMC, a Monte 
Carlo Markov Chain t ool, using the Chariot & Bruzual 2007 stel¬ 
lar population models (lAcauaviva et al.ll201 Ih . testing the effects 
of exponential star formation histories and more than one stellar 
population. 

We find 25 candidates which show a wide range of ages, 
masses and dust reddening. Our selection is designed to catch 
more massive and lumin ous galaxies than s i milar recent exam¬ 
ples in the literature ( e .g.lHagen et ^ 20141:IVargas et al.ll2014J: 
IWhitaker et al.l 120121: iGuaita et al.l l20Hlk which is expected 
given the relatively shallow depth but large volume probed. 

We see a tendency for the least massive systems to be dom¬ 
inated by younger stars, and vice versa (see Eig. [Tol l. This is 
usually seen in SED modelling when assuming a constant star 
formation rate (as in model A), since older gala xies are forced 
to be more massive (see e.g. [Vargas et al.l 1201411 . However, our 
models B and C indicate that most, if not all, of our LAE candi¬ 
dates show signs of recent elevated star formation, and that the 
light from the least massive systems is completely dominated by 
younger stars. 

Three out of nineteen (within the footprint) LAE candidates 
are found in the GOODS-Herschel catalog, showing far-infrared 
fluxes L/« > IO'^Lq. placing them in the LIRG-ULIRG regime. 
An estimate of the energy absorbed in the optical-UV shows 
these three candidates to have among the highest Luv,abs of the 
sample (Eig. as expected. We argue that these LAEs are ex¬ 
pected to be bright far-IR emitters based on their Luy^ats alone, 
and several of our remaining LAEs may be just below the detec¬ 
tion threshold in Herschel/PACS based on this energy budget. In 
conclusion, a lack of dust does not appear to be a requirement 
for Ly a escape. 
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Table 4. Best fit SED parameters, model A 


Candidate 

Quality 

Age (Gyr) 

low 

high 

E(B-V) 

low 

high 

Mass (10"^Mq) 

low 

high 


d.p. 

R-1 

NB392_6-1 
NB392 6-2 

N/N/I 

S/I/I 

0 16+**’^ 

0 

0.14 

0.05 

0.27 

0.12 

< 0.45 

0 74+0.08 
-0.13 

0.00 

0.35 

0.20 

0.38 

7 7+3.3 

■+-+_1.4 

14 6^^-^ 

3.1 

9.9 

3.7 

13.6 

0.18 

0.75 

12 

12 

0.004 

0.008 


NB401_6-1 

NB401_6-2 

NB401_6-3 

NB401_6-4 

NB401_6-5 

NB401_6-6 

NB401_6-7 

NB401_6-8 

NB401_6-9 

NB401_6-10 

NB401_6-11 

NB401_6-12 

NB401_6-13 

NB401_6-14 

NB401_6-15 

NB401_6-16 

NB401 6-17 


N/B/B 

N/N/I 

N/I/I 

N/I/I 

N/B/B 

S/I/I 

N/N/I 

N/B/B 

N/I/I 

N/B/I 

N/N/I 

N/N/B 

N/B/B 

N/I/I 

S/I/I 

N/N/B 

N/B/B 


0-90!n64 

n 70 + 0.14 

2.56™'| 

0 

^■^^-0.00 

A 99+0.46 
^ ~Q.QS 

O-lO-ni 
0-08!|” 
0-08!|| 
0 .10™-® 
0.84 


0.15 

1.61 


mi 


0.60 

0.04 

0.40 

0.05 

0.58 

0.72 

0.00 

0.17 

0.16 

0.22 

0.05 

0.04 

0.02 

0.07 

0.50 

0.07 

0.55 


2.99 

0.12 

2.99 

0.13 

2.99 

2.99 

2.99 

0.48 

0.29 

0.33 

0.15 

0.13 

0.07 

0.10 

2.99 

0.17 

2.99 


< 0.28 
< 0.28 

< 0.36 

< 0.33 

< 0.35 

< 0.23 

< 0.20 

< 0.33 

0 25+°°^ 

/tn 

0.21!; 
0.27!; 
0.21!; 
0.39!;, 

< 0.33 

< 0.26 
< 0.26 


04 


0.00 

0.00 

0.00 

0.00 

0.00 

0.21 

0.00 

0.00 

0.00 

0.21 

0.19 

0.24 

0.19 

0.38 

0.00 

0.00 

0.00 


0.12 

0.12 

0.18 

0.16 

0.15 

0.26 

0.11 

0.11 

0.14 

0.24 

0.21 

0.27 

0.22 

0.40 

0.14 

0.11 

0.12 


15 9+'5.4 

I 6+^-3 

II 

4 4 + 1.5 

9 . 3 !^-^ 

19 4+5.2 

015711:^3 

9 9+2.4 

^•^- 0.8 

14.6!« 

1 9 9+26.6 
6 9+'®'4 
1 4 + 0.9 

-1.3 

74 9+4.8 

7 . 9 !^-^ 

1 . 5 ^ 

6.4^! 


13.4 

0.8 

5.3 

2.7 
8.1 

10.7 
0.19 
1.1 

2.7 

14.3 
7.6 

4.1 
0.6 

20.1 
6.4 
0.9 
4.0 


18.1 

1.8 

7.7 

4.2 

10.3 
15.0 
1.05 

1.8 

3.2 

16.4 
15.6 

9.1 

1.2 
23.9 

8.7 
1.5 

5.7 


0.33 

2.15 
0.62 
2.03 

1.16 
1.39 
0.76 
0.22 
0.66 
0.37 
1.96 
1.18 
0.32 
6.21 
1.13 
1.52 
0.54 


10 

6 

11 

10 

12 

12 

6 

8 

11 

12 

6 

6 

8 

12 

9 
8 

10 


0.008 
0.002 
< 0.001 
0.048 
0.002 
0.002 
0.006 
0.005 
0.128 
0.003 
0.018 
0.009 
0.086 
0.002 
0.002 
0.004 
0.002 


NB401_ 

NB401_ 

NB401_ 

NB401 

_6-18 

_6-19 

_6-20 

6-21 

N/N/I 

S/B/B 

N/B/B 

N/N/B 

0.38!“^ 

0 17+ff23 

0.04 

0.16 

0.04 

0.02 

2.99 

0.22 

0.10 

0.14 

0 

^*■^^-0.14 

^•■^^-0.07 

A 99+0.09 
^•■^^-0.10 
<0.39 

0.29 

0.33 

0.33 

0.00 

0.34 

0.37 

0.36 

0.16 

22.6!iJ 

40.1 is 

5 4+r^-' 

-4-4 

1 6+®'2 
-14 

5.7 

39.2 

3.0 

0.8 

20.6 

43.3 

4.7 

2.6 

1.46 

0.54 

0.82 

1.30 

6 

12 

11 

6 

0.007 
0.003 
< 0.001 
0.056 

NB406 

6-l 

N/I/I 


0.32 

0.58 

< 0.19 

0.00 

0.08 

9 9 + 1.9 

^•■^-1.0 

1.9 

2.5 

1.29 

11 

0.004 

NB412 

18-1 

S/I/I 

0 19+456 

0.19 

0.31 

0 29+‘'‘)7 

0.24 

0.28 

-4.2 

20.9 

23.7 

0.30 

10 

0.015 


Notes. Derived SED properties for mo del A. For a detaile d description on the setup of the SED fit, see the text. The Quality column lists: Spectroscopically matching LAEs at the correct redshift 
according to the lBarger et alj ( l2008t) or l Adams et alj(l201 ih spectroscopic catalogs (S or N for Spectrum or No spectrum) / Possible confusion from HST imaging (B or I or N for Blended or Isolated 
or No data) / Possible confusion from ground-based imaging (B or I fo r Blend ed or Isolated). Mass here corresponds to the total mass that has gone into stars (i.e. the integral of the star formation 
rate over the age of galaxy). The R-1 value shows the lGelman & RuhinI ( Il992h R statistic. The best fit values and their upper and lower uncertainties correspond to the least-squares values with 68 % 
confidence intervals. If only an upper or lower limit is listed, this corresponds to the 68 % limit. The two columns following each parameter give the 68 % lower and upper limits when marginalized 
over all parameters. represents the 4 ^“ of the best fit model, divided by the degrees of freedom, which is the number of data points (d.p.) minus the number of parameters fitted. 
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Table 5. Best fit SED parameters, model B 


Candidate 

Quality 

E(B-V) 

low 

high 

Mass (lO^Mo) 

low 

high 

T -' (Gyr-') 

low 

high 


d.p. 

R-1 

NB392. 

_6-l 

N/N/I 

< 0.36 

0.00 

0.15 

9 0 + 4.1 
^•^- 0.9 

3.4 

4.2 

< -0.30 

-6.34 

-1.99 

0.30 

12 

0.015 

NB392. 

6-2 

S/I/I 


0.25 

0.28 

13 9+^' 

13.4 

15.5 

< -1.70 

-6.34 

-2.77 

0.97 

12 

< 0.001 

NB401. 

_6-l 

N/B/B 

< 0.23 

0.00 

0.12 

21 

15.4 

19.5 

-0 30+" '" 

u .^ u _2 45 

-1.45 

-0.60 

0.43 

10 

0.097 

NB401. 

_6-2 

N/N/I 

< 0.22 

0.00 

0.08 

2 1 +^® 
lll +^-3 

46+^^" 

112+^-3 

2.5 

3.4 

< - 1.10 

-6.34 

-2.52 

2.45 

6 

0.001 

NB401. 

_6-3 

N/I/I 

< 0.28 

0.00 

0.22 

5.6 

7.6 

< 1.03 

-4.94 

-0.36 

0.56 

11 

< 0.001 

NB401. 

_6-4 

N/I/I 

< 0.22 

0.00 

0.08 

4.2 

5.0 

< -2.27 

-6.35 

-3.25 

2.12 

10 

< 0.001 

NB401. 

_6-5 

N/B/B 

0.13!“;® 

0.13 

0.17 

8.6 

10.8 

-0.94!'-40 

0 33+'-® 

+'■-’■+- 4.83 
-0 49 + 0-23 
+'■^'^-5 83 

-2.87 

-0.66 

1.26 

12 

0.011 

NB401. 

NB401. 

_6-6 

_6-7 

S/I/I 

N/N/I 

< 0.51 

< 0.17 

0.00 

0.00 

0.28 

0.10 

19 2+2-^ 

74-^" 

'*^60 

12.0 

2.9 

18.1 

4.4 

-2.92 

-6.32 

0.86 

-1.42 

1.38 

1.20 

12 

6 

0.003 
< 0.001 

NB401. 

_6-8 

N/B/B 

< 0.18 

0.00 

0.10 

2 5+^-^ 

1.4 

1.9 

< -0.27 

-5.45 

-1.67 

0.23 

8 

0.008 

NB401. 

_6-9 

N/I/I 

< 0.29 

0.00 

0.12 

9 q + 2.6 

3.1 

3.7 

< -0.83 

-6.34 

-2.48 

0.65 

11 

0.002 

NB401. 

_6-10 

N/B/I 

0 94 + 0.04 

0 19+0® 

0 08 

0.23 

0.25 

15.4 

18.3 

24 

-6.34 

-2.59 

0.37 

12 

< 0.001 

O O 

_6-ll 

_6-12 

N/N/I 

N/N/B 

0.17 

0.20 

0.19 

0.23 

15 6+^^-^ 

I'+t 

20.1 

13.9 

26.6 

18.6 

< -1.28 

-6 13+^-^’^ 

^ 0 00 

-6.33 

-4.50 

-2.46 

-1.32 

2.17 

1.55 

6 

6 

< 0.001 
< 0.001 

NB401. 

_6-13 

N/B/B 

< 0.32 

0.00 

0.13 

1 3+^-^ 

04 

1.5 

1.9 

< -1.38 

-5.83 

-2.84 

0.65 

8 

< 0.001 

NB401. 

_6-14 

N/I/I 

0 32+0 02 
-0-03 

0 12+0'' 

0.30 

0.32 

18.6+1',2 

21.1 

24.0 

< -4.61 

-6.35 

-4.94 

7.65 

12 

0.023 

NB401. 

_6-15 

S/I/I 

0.13 

0.18 


6.5 

8.4 

_l 19 + 1.83 
+ ■ + ^- 4.94 

-4.36 

-0.97 

1.16 

9 

0.009 

NB401. 

_6-16 

N/N/B 

< 0.19 

0.00 

0.07 

1 4 + 1.5 
^‘^05 

1.5 

2.0 

< -1.40 

-5.82 

-3.06 

1.56 

8 

0.002 

NB401. 

_6-17 

N/B/B 

< 0.27 

0.00 

0.16 

6 7+^' 

+’•'- 4.0 

99 4 + 46.6 

38:6^5 

4 5+3-6' 

4.0 

5.3 

-0.27:°-™ 

-4.26 

-0.71 

0.54 

10 

< 0.001 

NB401. 

NB401. 

_6-18 

_6-19 

N/N/I 

S/B/B 

0 

^•'^^- 0.20 

0 36+° °^ 
^•*^^-0.06 

QQJ 

0.27 

0.33 

0.32 

0.35 

18.7 

41.7 

27.7 

48.3 

<0.87 
< -2.50 

-5.29 

-6.34 

-1.02 

-3.48 

1.46 

0.56 

6 

12 

0.005 
< 0.001 

NB401. 

_6-20 

N/B/B 

0.24 

0.27 

5.0 

5.9 

< -2.32 

-6.34 

-3.21 

0.94 

11 

< 0.001 

NB401. 

6-21 

N/N/B 

< 0.31 

0.00 

0.13 

3.5 

5.0 

< -0.33 

-5.57 

-1.86 

1.40 

6 

< 0.001 


NB406 

6 -I 

N/I/I 

< 0.19 

0.00 

0.09 

9 7 + 2.0 

^•'- 1.5 

2.0 


-4.97 

-1.41 

1.34 

11 

0.013 

NB412 

18-1 

S/I/I 

0 29+''-'" 

_+'■^+^- 0.09 

0.25 

0.28 

21 5+'++’ 

22.4 

26.4 < -1.43 

-5.72 

-2.61 

0.32 

10 

0.003 


Notes. Same as Tabled but for model B. The Age of each candidate is fixed at 2.2 Gyr. 
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Table 6 . Best fit SED parameters, model C 
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10 
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10 
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Notes. Same as TablelH but for model C. fymmg represents the fraction of the mass found in the younger population. A missing value means the fraction could not be constrained. 
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Appendix A: Postage stamps of all candidate LAEs 

We here show the postage stamps for all of our candidates in 
Figure IaI 


A&A-aadem, Online Material p 18 



Fig. A.l. Postage stamps for all of our candidate LAEs. The field of view of all stamps is 10 x 10 arcsec. The large RGB picture shows an HST bvi 
color composite made with the scale i n each color channel related to the true total intensity, with an arcsinh intensity scaling relation. This coloring 
scheme is similar to that described in lLunton et alJ ll2004h . The white circle has a one arcsec radius - typically twice the size of the ground-based 
seeing disk. 
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Fig. A.2. Stamps (cont.) 
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Fig. A.3. Stamps (cont.) 
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Fig. A.4. Stamps (cont.) 
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Fig. A.5. Stamps (cont.) 
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Filter 401_6, candidate 17 
12:36:59.57 +62:19:37.7 
U-NB: 1.14 ± 0.93 
g'-NB: 0.61 ± 0.94 
B-Z: -0.06 ± 0.40 
Z-Ks: 1.03 ± 0.48 


Filter 401_6, candidate 20 
12:37:24.47 +62:18:18.6 
U-NB: 1.58 ± 0.70 
g'-NB: 0.35 ± 0.68 
B-Z: 0.74 ± 0.33 
Z-Ks: 1.42 ± 0.33 


5.8 Tim 


8.0 iim 




Filter 401_6, candidate 18 
12:36:26.20 +62:19:32.1 
U-NB: 1.28 ± 0.99 
g'-NB: 0.21 ± 0.99 
B-Z: 0.57 ± 0.23 
Z-Ks: 1.71 ± 0.28 


Filter 401_6, candidate 19 
12:37:20.10 +62:19:23.6 
U-NB: 1.31 ± 0.58 
g'-NB: 0.60 ± 0.58 
B-Z: 0.87 ± 0.13 
Z-Ks: 1.78 ± 0.22 
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Fig. A.6. Stamps (cont.) 
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Filter 401_6, candidate 21 
12:38:21.32 +62:17:50.4 
U-NB: 1.96 ± 0.81 
g'-NB: 2.87 ± 1.07 
B-Z: 0.35 ± 0.23 
Z-Ks: 0.36 ± 0.32 


Filter 406_6, candidate 1 
12:37:11.29 +62:17:09.5 
U-NB: 2.44 ± 0.80 
g'-NB: 1.91 ± 1.55 
B-Z: 0.23 ± 0.27 
Z-Ks: 0.02 ± 0.37 
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Filter 412_18, candidate 1 
12:36:32.05 +62:12:52.7 
U-NB: 1.42 ± 0.96 
g'-NB: 0.60 ± 0.96 
B-Z: 0.63 ± 0.15 
Z-Ks: 1.45 ± 0.24 


5.8 pm 


Fig. A.7. Stamps (cont.) 
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